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TECHNICAL MEMORANDUM 


DISTRIBUTED CONTROL USING LINEAR MOMENTUM 
EXCHANGE DEVICES 

I. INTRODUCTION 


This report documents Marshall Space Flight Center's (MSFC) role in the Air 
Force Wright Aeronautical Laboratory's (AFWAL) Vibration Control of Space Structures 
(VCOSS) experiment. This experiment, often referred to as VCOSS-II, was designed 
as a hardware demonstration stemming from the original VCOSS program which, in 
turn, was an outgrowth of AFWAL's Active Control of Space Structures (ACOSS) pro- 
gram. The purpose of ACOSS was to develop control strategies for Large Space 
Structures (LSS); VCOSS investigated the means to implement the ACOSS algorithms, 
while VCOSS-II concentrated on a hardware demonstration of VCOSS control elements. 
In particular, VCOSS-II was the active vibration suppression of an LSS test article 
by means of Linear Momentum Exchange Devices (LMEDs) . Several LMEDs were to be 
strategically located on a generic LSS test structure to sense any vibratory motion 
of the structure and transfer the associated linear momentum to a proof mass using 
active control elements. The LMED controller, in turn, would dissipate the motion 
of the proof mass, thereby increasing the overall damping of the LSS test article. 

Since NASA and the Air Force have mutual interests in the active control of 
LSS, the MSFC Ground Facility for Large Space Structures Control Verification 
(GF/LSSCV) was selected as the experimental site for VCOSS-II. AFWAL awarded a 
contract to the TRW Space and Technology Group to design, fabricate, deliver and 
demonstrate LMED hardware on the GF/LSSCV. A thirteen meter ASTROMAST with 
an 8- in. triangular cross-section was selected as the test article. This lightweight, 
flexible beam, which was formerly a spare magnetometer boom for the Voyager space- 
craft, has previously been used at MSFC for LSS control experiments. The beam is 
attached vertically to a three degree of freedom, flight quality pointing mount, with 
two complete Inertial Reference Units (IRUs) attached one to either end of the beam. 
The LMEDs were designed to complement this assortment of control hardware. 

Planning included hardware demonstrations of the LMEDs alone and in conjunction 
with the GF/LSSCV three gimbal pointing mount. 

The original VCOSS-II test plan was arranged along the following lines. MSFC 
agreed to provide the test facility, perform preliminary modal tests on the test 
article and provide mechanical integration support for the experiments. MSFC agreed 
also to provide an independent preliminary structural dynamic model verified by the 
modal tests. In exchange for the use of the GF/LSSCV, AFWAL granted MSFC 
unrestricted use of the VCOSS-II hardware for an indefinite period to perform 
independent control studies, provided that these studies did not interfere with the 
VCOSS program schedule. 

This report will address the MSFC control experiments. A brief description of 
the facility, control hardware and test configuration is followed by a summary of 
dynamic models of the structure and of the LMEDs. This will be followed by a 
description of the MSFC control strategies. Experimental results are then presented. 

A brief section of conclusions and recommendations appears at the end of the report. 

It should be stressed that this report is intended to document the MSFC work to date, 


and not to provide a comprehensive analysis of the applicability of LMEDs for LSS. 
The work to be presented reflects a status report of the initial testing performed 
through early 1986. Experimental work with the LMEDs is continuing, the results of 
which will be published in future reports. 


II. EXPERIMENT DESCRIPTION 


This section will describe the experimental facility, the hardware components 
and the experimental configurations. Dynamic models of each configuration developed 
by the Control Dynamics Company for MSFC will be discussed. Analytic models of 
the basic LMED components will also be discussed. 


2.1 Experimental Facilities 

2.1.1 GF/LSSCV 

Figure 1 outlines the basic components of the LSS test article in the GF/LSSCV. 
The ASTROMAST is shown suspended vertically from a modified version of the Sperry 
Corporation's engineering model of the Advanced Gimbal System (AGS). The AGS 
consists of two gimballed dc torque motors capable of developing 51 Newton meters 
(N-m). The AGS is mounted below an air bearing about the vertical axis with a 
19 N-m torque motor. All three torque motors have measured bandwidths in excess 
of 100 Hz. The modified AGS is attached to a hydraulically driven, two translational 
degree of freedom, Base Excitation Table (BET). A three-axis accelerometer unit 
is located on the BET, and an identical accelerometer unit is located on the tip sensor 
package (at ground level). These accelerometers have digital outputs, 25 Hz band- 
width and 11 micro- g sensitivity. Three single axis rate gyros are attached to the 
AGS mounting plate. These gyros have 40 Hz bandwidths, 2 arc- second /second 
sensitivity, one degree per second saturation limits, and provide analog output 
voltages. The tip rate gyros have 70 Hz bandwidths, digital outputs, and 50 to 90 
arc- second /second sensitivity. A MSFC built COSMEC computer performs the 
Input /Output (I/O) task of interfacing the sensors and actuators with a HP- 9000 
control computer. The HP- 9000 and COSMEC are currently configured to operate at 
a 50 Hz sample rate. The HP- 9000 performs digital sensor indication updates, coor- 
dinate transformations, IRU functions and control law implementation, all of which is 
accomplished in the HP BASIC language. The ASTROMAST is 13 m long with a 20.3 
cm on a side equilateral triangular cross-section and three continuous S- glass 
longerons. There are 92 flexible battens distributed horizontally along the beam in 
14.3 cm increments. These locations will be referenced as station numbers with 
station 1 at the AGS mounting plate and station 92 at the tip. The ASTROMAST has 
a static twist of 260 deg about the vertical axis when loaded as in this configuration. 
Figure 2 is a photograph of the Astromast in a previous experimental configuration. 

A 4-m cruciform is shown attached to the tip sensor package. 


2.1.2 Linear Momentum Exchange Device 

Two pairs of LMEDs are identified in representative locations in Figure 1. 

Figure 3 is a close-up photograph of one such LMED pair mounted to the ASTROMAST 
The schematic drawing in Figure 4 identifies the major components. Two LMEDs are 
shown mounted orthogonally to each other on opposite sides of a rigid aluminum frame 
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Figure 1. MSFC LSS Ground Test Verification Facility 







Figure 2. Photograph of GF/LSSCV ASTROMAST. 
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Figure 3. Close-up photograph of LMED mounted on ASTROMAST. 
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Figure 4. Schematic drawing of LMED construction. 


which substitutes for one of the ASTROMAST flexible battens. Each LMED consists 
of a 0.75 Kg cylindrical moving mass with ±1.27 cm of linear travel along a stainless 
steel shaft; a dc coupled, electronically driven voice coil attached to the frame; a 
Linear Variable Displacement Transducer (LVDT) which measures the position of the 
proof-mass relative to the frame; and a linear accelerometer with a pre-amplifier 
stage. The input axis of the accelerometer is aligned with the stainless steel shaft. 
The proof- mass contains samarium cobalt permanent magnets which move outside of 
the actuating coils, and stainless steel ball bearings which support the moving mass 
internally. The LVDT core is teflon coated, and its attachment to the moving mass 
restricts rotational motion. The moving masses of the two LMEDs are placed at the 
centroid of the triangular frame. The toteil non-moving mass of the LMED pair is 
1.0 Kg which gives a total mass of 2.5 Kg for a complete pair of devices. As an 
option, each LMED can employ centering springs on either side of the proof-mass. 
Figure 3 also shows laser diode sources mounted on ’’diving board” extensions. 
These lasers are part of an optical position sensing system which was not used in 
VCOSS-II, and will not be discussed in this report. 


2.2 Experimental Configuration 

The VCOSS experiment ultimately involved three test configurations as shown 
schematically in Figure 5. In the first configuration identified as configuration 1-A, 
the AGS and BET were rigidly constrained and the 18-Kg tip sensor package was 
removed. The two LMED pairs were located at stations 91 and 49. These locations 
were selected for the following reasons. The middle LMEDs, which were labelled A 
and B, would have control authority throughout the first 20 bending modes, especially 
at the mid-frequency modes. The lower pair of actuators would have control 
authority over the first eleven bending modes, primarily at lower frequencies. In 
addition, the 42 station separation corresponds to 120 deg of twist of the ASTROMAST. 
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Figure 5. Schematic of the VCOSS-II test configurations 


Since the LMED frame is an equilateral triangle, the middle device could be reoriented 
by 120 deg to align it with the lower device. The discrete possibilities for mounting 
these devices implies a 60 deg (21 station) separation intended for best alignments; 

30 deg intervals will introduce ±1.4 deg offsets. A spare LMED (a misnomer, since 
this device lacks an accelerometer) was located at station 22 to serve as an excitation 
source. This device, labelled LMED-E, was oriented about 43 deg relative to the 
other actuators. To prevent the relatively massive actuators from creating mid- 
frequency vibration nodes on the structure, three 1.55 Kg lumped masses were 
attached to the structure (one each longeron) at station 36. 

Configuration 1-A, described above, allowed the LMEDs to be tested independent 
of the AGS system, thereby simplifying the modeling task. Configuration 1-B is 
similar to 1-A except that the AGS system was free to rotate, which caused the fre- 
quencies of vibration to decrease. In configuration 2, the tip sensor package was 
mounted to the ASTROMAST along with a 3.6 Kg aluminum cruciform, with four 
approximately 2-m long, constant cross-section (0.635 cm) arms. This had the com- 
bined effect of further decreasing the modal frequencies while increasing the modal 
density, and coupling the bending and torsion modes. The excitation source, 

LMED-E, was removed along with the three lumped masses. The middle pair of 
actuators, LMEDs A and B, were located at station 13 and the lower pair at station 
55. The BET was unlocked and could serve as the excitation source, as could the 
AGS torquers. 

The three configurations represent a departure from the original test plan. 
Configuration 1-A was planned as the baseline to demonstrate and characterize the 
LMED hardware components. The addition of the aluminum cruciform to the configura- 
tion was planned to demonstrate the robustness of the controller. Finally, the AGS 
gimbals were to be unlocked and actively driven along with the LMEDs to demonstrate 
multiple input /multiple output control techniques. As it turned out, configurations 
1-A and 1-B were used to test an analog controller. Configuration 1-B was also used 
to demonstrate the MSFC designed digital control of the LMEDs. Configuration 2 was 
used for a demonstration of the simultaneous use of the AGS torque motors and the 
LMED systems for control. This selection of test configurations was dictated by 
schedule constraints which, combined with GF/LSSCV hardware problems, required 
unrestricted access to the system components for troubleshooting and testing. 


2.3 Structural Dynamic Models 

MSFC and the Control Dynamics Company spent a considerable amount of time 
generating structural dynamic models of the planned VCOSS test configurations and 
verifying these models with extensive modal tests on the structure. However, the 
last minute changes in these configurations (as described in the previous section) 
meant that the corresponding changes to the dynamic models could not be verified by 
modal tests. Thus, the models to be presented herein provide reasonable levels of 
accuracy, but not necessarily the best that could be achieved given full experimental 
information. 

Figure 6 depicts the elements used to model test configuration 1-A. The local 
coordinate frames associated with each element define the numerical degrees of freedom 
of the overall structure. The first index along an axis indicates a translational 
degree of freedom; the second index indicates a rotational degree of freedom. The 
analysis was performed with the ISMIS (Interactive Structures and Matrix Interpre- 
tive System) finite element code. This model contains 41 degrees of freedom including 
the five moving masses associated with the five LMEDs. The LMED models also 
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included centering springs which produced a set of five nested modes near 1 Hz. 
Table 1 provides the frequencies and the mode shapes at each device. For example, 
LMED-E exerts a force on the moving mass at degree of freedom number 37; the 
components of an equal and opposite force are applied to the structure at locations 7 
and 8. Torsional modes in this model, which is highly symmetric, are evident in 
Table 1 by the mode with zero lateral displacements, i.e., modes 10, 13, 18, 19, 20, 
23, and 24. These modes can be neglected in subsequent analyses since they do not 
interact with the LMEDs in this model. The last column on the right identifies the 
frequencies of vibration which were evident in the open loop frequency response data 
presented in Section IV. Figure 7 and Table 2 provide the analogous information for 
configuration 1-B; Figure 8 and Table 3 apply to configuration 2. 
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Figure 6. Dynamic model of configuration 1-A. 


9 



TABLE 1. VCOSS-II DYNAMIC MODEL FOR CONFIGURATION 1-A 
(AGS Gimbals Locked; No Tip Sensors or Cruciform) 
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Figure 7. Dynamic model of configuration 1-B. 
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TABLE 2. VCOSS-II DYNAMIC MODEL FOR CONFIGURATION 1-B 
(Gimbals Unlocked; No Tip Sensors or Cruciform) 
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Figure 8. Dynamic model of configuration 2. 



TABLE 3. VCOSS-II DYNAMIC MODEL FOR CONFIGURATION 
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2.4 LMED Analysis and Modeling 

This section will describe a limited dynamic model of the LMED including analy- 
sis of the linear motor drive electronics, actuator dynamics, the accelerometer and 
accelerometer preamplifier, and the LVDT relative position sensor. The model is 
developed analytically and supported with laboratory measurements. A high fidelity 
model suitable for computer simulation was not completed due to the complexity and 
magnitude of the measured non-linearities. The level of influence of these non- 
linearities over the actuator dynamics was later greatly diminished by techniques 
developed in the MSFC laboratories in conjunction with Control Dynamics Company. 

Figure 9 provides a transfer function description of the major LMED components. 
This model represents a single device as tested on a laboratory bench, for example, 
on an ideal air bearing, i.e., "slippery plate." The forward path describes the 
input /output relationship with respect to the experiment control computer; V cm( i is 
the command voltage from the controller and V a is the return signal from the accel- 
erometer preamplifier. The feedback path describes the dynamics of the motion of 
the "proof mass"Mp, incorporating viscous damping and centering springs. The sum 
of these feedback return forces with that produced by the motor's solenoid coil is 
applied in an equal but opposite manner to the proof mass and the actuator base. 

The following parameters are identified in the model. 


K D = 0.2 (A/V) the frequency independent gain of the driver electronics 
= 2 x 10 4 (sec) 

= 4 x 10 ® (sec) 

K m = 2.1 (lb/A) the dc motor gain 

= 1.3 (mA/g) the dc gain of the accelerometer 

C a = 0.3 to 0.8 specified range of damping ratio of accelerometer 

w = 2 tt 500 (r/s) specified minimum natural frequency of the accelerometer 
8. 

Kp = 30 (V/A) dc gain on accelerometer preamp 
t 3 = 0.9 (msec) 

Mp = 1.65 (lbm) mass of the moving mass 

Mg = 3.85 (lbm) mass of base and non-moving (orthogonal) proof mass 


15 



Figure 9. Block diagram of LMED analytical model. 


D = 0.046 (lbf-s/in.) estimated damping coefficient 


K =0.2 (lbf/in.) the combined linear centering spring coefficient 
s 

K r = 20 (V/in.) the LVDT gain 

V cm d = -10 (V) the range of input command voltages 
i = (A) current into motor coil 

F c = (lb) force produced by motor coil 

i = (A) output current from accelerometer 
a 

V = ±10 (V) range of accelerometer preamp return signal 
a 

V R = ±10 (V) range of output voltage from LVDT 





Xp = (in.) motion of the proof mass in laboratory frame 

Xg = (in.) motion of the base in laboratory frame 

X 0 = (in.) motion of proof mass relative to base . 

K 


Notice the motor dynamics do not include the effects of back emf produced by 
the velocity of the proof- mass permanent magnets relative to the motor coils. This 
effect is negligible due to the constant current feedback amplifier which drives the 
motor. Also note that the mass of the base includes the proof mass of the actuator 
mounted perpendicular to the modeled device. 

There are no dynamics associated with the LVDT in the model since it was not 
used in the control loop. The LVDT data sheet specifies a 200 Hz cutoff frequency 
with less than 0.25 mV of ripple and 0.25 percent linearity. 

The drive electronics transfer function was obtained from the schematic shown 
in Figure 10. Assuming an ideal operational amplifier (infinite gain high input impe- 
dance and low output impedance) the following transfer function between the common 
voltage and motor current can be obtained. 


m 


A o <R L + Z F* " R L 


V cmd A o *4. R I + [<Z m + «I + R L + V + R L (R I + V 1 


where A q = ideal amplifier gain > 10 


R t 


Z F 1 + C F R p s' 


Z = 


( r M1 + 1/C M s) (R M£ + L M s) 


m (R M1 + R M2> + L M S + 1/C 1 S 


Y k b k m s 


M b - Mp 


M p Mg s' + (M p + M fi ) Ds + (M p + M fi ) K g 


Kg = 0.237V /in. /sec back emf constant 

The transfer function I m / V cmd represents a constant current amplifier. Since A Q is 
large , this relation can be closely described by the approximation 
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^cmd 


R^+ Zp _ -K d (1 + x^s) 
R L R I = 1 + T 2 S 


The maximum error in this approximation is 2 percent which occurs at the 40 kHz 
resonance of Z . The term y is highly damped so it has little influence compared to 

A q . It does, however, encompass the back emf effect which shows that back emf 

need not be included when using a constant current source with high gain and high 

input impedance. Furthermore, since ~ 4 x 10 sec, the drive electronics can be 

sufficiently represented by K D at the structural frequencies, and the conversion from 

command volts to force becomes K m = 0.42 (lbf/V). 

The accelerometer preamplifier circuit is shown in Figure 11. The ideal opera- 
tional amplifier assumption yields the transfer function 



1 

1 + t 3 s 


where = 0.9 msec. The 177 Hz break frequency of the amplifier con- 

tributes only eight degrees of phase lag at 25 Hz, and may sometimes be neglected. 
Similarly, the 500 Hz low pass effects of the accelerometer may be neglected in some 
cases. 



Rj=200kft 


Figure 10. 


R F =40Kft R m =l.Q C M =0.01yf R L =1ft 

Cp=0.001pf R ri2 =6ft L M =1.52mh 

Kg=0.237v/in/sec 

Schematic of LMED drive electronics. 
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R 1 



(to controller) 


R 0 =lkfl R-|=30kfi R2=100fl 


Figure 11. Accelerometer output amplifier schematic. 

Figure 12b provides measured frequency response data for LMED C, taken in 
the laboratory. The LVDT response was obtained with the LMED base clamped to the 
laboratory test bench. The accelerometer response was obtained with the LMED float- 
ing on an air bearing. Figure 12a shows the predicted response using the ideal 
analytical models. The LVDT response compares well with measured data; the accel- 
erometer response does not match measured data very well. A subsequent report 
from Control Dynamics to MSFC, entitled ’’Characterization and Hardware Modifications 
of LMEDs," describes the unmodeled nonlinearities of the LMEDs and the steps taken 
to linearize them. The hysteresis and stiction initially inherent to the LMEDs lead to 
the discrepancies seen. 


III. CONTROLLER DESIGNS 


3.1 MSFC Digital LMED Controller 

To gain familiarity with the performance capabilities of the LMEDs, MSFC adopted 
a conservative control methodology. Since the vibration modes in Experiment 1 were 
clearly separable in frequency, a narrow bandwidth, digital, rate feedback controller 
was designed for each LMED. The center frequency of each controller was placed 
near the frequency of the dominant mode at that actuator location. Referencing back 
to Table 2, LMED A was tuned to the 1.96 Hz mode; LMED B to 2.78 Hz, LMED C to 
11.8 Hz and LMED D to 12.88 Hz. Figure 13 provides the basic block diagram of the 
control algorithm. The HP-9000/COSMEC control computer samples the accelerometer 
pre-amplifier signals at a 50 Hz rate; the A/D input filters on the COSMEC were 
bypassed. The second-order digital filter performs both discrete time integration and 
bandpass filtering of the accelerometer signal. The narrow bandpass filter helps 
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Figure 12a. Analytic frequency response from LMED to accelerometer 










Rate Comp. 

Integrator High Q Bandpass Gain Delay 



COSMEC 

Figure 13. Block diagram of LMED digital controller. 

reduce the effects of the LMED non-linearities. Note that in Figure 13 a unit delay 
z 1 accounts for a fixed computational delay of one sample period. 

The filter coefficients are determined as follows. A second-order Butterworth 
bandpass filter is described by 

(w /Q)s (w„ - u,)s 

G B p (s) = ~2 2 2 = S L 

s + + u o s + ( w 2 ~ “f) S + to 1 0) 2 


where 

to = 2nf o = center frequency 

2 

Q = u) 0 /Aa) Uj = (o Q /a>2 


Aw =u) 2 ~u)^ = 3db bandwidth 


- w o (1 


+ A 


1 + 4 Q )/(2Q) 


This continuous filter is transformed to discrete time using the Tustin transformation 
with frequency prewarping; 


_ 2 z-1 . i _ 2 , .1 . 2 , 

s ^ ^ tan ( 2* Uj^T) ^ tan 

where T is the sample period. Substituting for s and replacing £1. , n 0 
results in 


4 " 2 T > 

for Wj, a> 2 
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where 


Ggp(z) “ 


b 1 ( z 2 ~ 1) 


a 0 z + a, z + a 
2 1 o 


b l = (y 2 " Y l> 
a 2 = 1 + Y 2 ' Y 1 + Y 1 Y 2 
a l = 2 (Y 1 y 2 ' X) 

a o = 1 + y i - y 2 + Y i y 2 

= tan u^T) = ^ ftjT 


Y 2 = tan (i go 2 T) = ^ ^ 2 T 


The discrete time integration is accomplished by replacing 1/s with 


T £+1 
2 z-1 


To implement the integrator filter with the control computer, the following reduction 
was made. 


U(z) 

Y(z) 


T z+1 
2 z-1 


b^z-lKz+l) 


a 2 z +a lZ +a Q 


b 2 + \ 2-1 + b 0 
a 2 + a l Z_1 + a o 





1 * 


Taking the inverse Z transformation yields the discrete time algorithm 


Ui=K r (b 2 Y i + bl + b Q Y,_2 « x 


a o u i-2 )/a 2 


where represents the desired rate feedback gain constant. The control algorithm 
was programmed so that K r , u> Q and Q for each LMED could be specified at runtime. 
The following values were eventually used during experiment IB . 
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f 

o 

Q 

K r 

LMED-A 

5 

10 

-750 

LMED-B 

5 

10 

-600 

LMED-C 

12.5 

10 

-750 

LMED-D 

12.5 

10 

-600 


3.2 MSFC Digital LMED/AGS Controller 

In Experiment 2 the LMEDs were operated in parallel, i.e., simultaneously with 
the AGS torque motors. Because this experiment configuration closely resembled a 
previous MSFC experimental setup, a previously developed, low bandwidth AGS con- 
troller was modified to work with the LMEDs [13. Although the resulting control law 
was functionally decentralized, the individual controllers were cast into state variable 
form to facilitate future analysis and investigation of other control methods. This 
particular control experiment was designed to demonstrate the functioning hardware 
and simple control techniques. The control gains and frequencies were selected 
largely experimentally to achieve a working closed loop system using seven control 
effectors and seven control sensors. The performance which was obtained was not 
fully supported analytically due mostly to the inherent non-linearities of the LMEDs 
which complicate such analyses. The following controller formulation, then is pro- 
vided for future control studies involving more linear actuators. A description is 
also provided of how the control gains and frequencies were selected experimentally. 

The discrete time filter given in the previous section can be expressed in the 
following state variable form 


[* X il 

cl 


0 -a 2 ' 


fx . I 

cl 


V 


= 




+ 


- X c2. 

■J _L 1 

- 1 '“I. 


CM 

o 

X 


. 6 2- 


where 

a 2 = V a 2 h = (b o " a o b 2 )/a 2 

«i = * 1^2 P 2 = ^ b l " a l b 2^ /a 2 

and Y C j = K r [ x c2 . + (b 2 /a 2 ) Y.] where X g . represents the accelerometer preamplifier 
signal, and Y c ^ is the control law output. 

The structure, sensor, actuator and controller can be represented in the follow- 
ing state variable form: 
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structure 

X^ = AX + BU 
P P P P P 

X = (2n xl) plant state vector 
P P 

measurement 

Y = C X + DU 
P P P P P 

Y = (mxl) plant output vector 

Jr 

sensors 

X c = AX + BY 
s s s s p 

X = (n xl) sensor state vector 
s s 

actuators 

U = A U + BY 
p up U c 

U = (n xl) plant input vector 

r U 

control gains 

Y = C X + D X 
c. C c. c s. 

Y = (£xl) controller output vector 
c i 

control filter 

X o = AX + B X 

c i+l c c i c s i 

X = (n xl) controller state vector 

c. c 


Note that the control law is written in discrete time coordinates, where X 0 is the ith 

s i 

sample of sensors X g . X is the controller state at the ith sample and is the 

controller output determined at the ith sample. The relationships between these state 
equations are shown in Figure 14. The plant equations are most easily described in 
modal coordinates, in which n represents the number of vibration modes in the plant 
model ; ” 


X p “ tn l V** ^np ^n^ 

where r\. represents the ith mode; A is block diagonal with each modal block A 

I P P. 

given by 



1 

-Vi . 


where go. and are the ith modal frequency and damping; and B is given by 

A Jr 


0 

b lll 

0 

b 12 1 

0 

bl "p 

0 

b 21 * 

0 

b 22 1 

0 

b 2n 

• 

• 

1 

l 


I 

| 

• 

P 

• 

0 

1 

b n' 

0 

1 

b . 2 l 

0 

b -p- 
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Figure 14. State space model of digital AGS/LMED controller and plant 
















where b.. = jth mode shape at the ith actuator. The matrices C and D are obtained 

J * nn lr 

as follows: Let Y = C, x +C„X , where C. and C„ each resemble B . In particular, 

plp2p lz 

for each mode j, j = l,2,...n , each (2j-l) column of contains mode shapes at 

each proportional measurement, each (2j) column of contains mode shapes at each 

rate sensor, and each (2j) column of C2 contains mode shapes for each acceleration 
measurement ; otherwise , the elements of C ^ and C 2 are zero . Substituting 

+ U into the above description for Y gives Y_ = C X + D U where C = 
PP P PPPPP P 

(C^ + C 2 Ap) and Dp = C 2 Bp. Taking advantage of the inertial reference data provided 

by the HP- 9000, provides angular position information from the base gyros while 

C 2 provides the accelerometer information from each LMED. The sensor dynamics 


provide first order approximations for the gyro (40 Hz bandwidth) and accelerometer 

(50 Hz) electronics. Therefore, n = m, and both A and B are (mxm) diagonal 

s s s 

matrices, with the negative of the sensor break frequencies (radians /second) along 

A and the sensor gain (volts /physical units) along B . The actuator dynamics need 
s s 

only represent the first order models of the AGS torque motors which have 100 Hz 

bandwidths. The LMED bandwidths exceed 2 KHz, which implies n y = i. For some 

analyses even the AGS dynamics can be neglected since the phase lag is less than 
12 degrees at 20 Hz. The discrete time control law is also block diagonal, with each 
block resembling the (2x2) control law given at the beginning of this section; this 
implies that n = 2m. 


For control law design and analysis it is convenient to discretize the continuous 
time equations. To do this, consider the following augmented state equation 



B C 
s p 


0 B 

A s 
0 


P 

B D 
s I 

A u - 




'x 


"o 



p 





X 

+ 

0 

) 


s 





U 


B . 

— 

1 u p 

u 


[Y c 3 


which can be written X = AX + B U. The system can be discretized for a sample 
period T, taking into account the constant computational delay of one sample period. 
The resulting system is given by 


X 


i+1 


A A 

= AX. + BU. 


where 
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and 


A 

B 


A 

-f 


e^ T diB 


Defining additional states V. = U. this can be expressed as 


" x i + r 


’ A 

B 


X. 

l 

+ 

0 

- v i + i . 


_ 0 

0_ 


V. 

1 J 


- 1 n - 

n 


U i 


The discrete time measurement equation becomes 


[X ] = [ C 0] 
S i 


X. 

.1 


V. J 


where c, = [0 2n I ns ° nu 3 » ! ns = identity matrix, and 0 2n , 0 Ru are zero matrices, 
Combining the control filter equations with the above gives 


X i + 1 


ABO 


X. 

l 

V i + 1 

- 

0 0 0 


V. 

l 

X 

c i+r 


.Be 0 A 

c c 


.X 

c 


0 

I 

] 

L0 


nu 


U i 


or X. +1 = A Xj + B U. 


Y. = CX. where C 


U 0 I J 

nc 


U. = KY. where K = [D C ] 
11 c c 


This places the entire system in standard discrete time, state variable form. 
The elements of the gain K and the filter element A c and B c can be modified at the 
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control designer's option. In Experiment II, these elements were determined by 
examination of measured open loop frequency response data. The open loop transfer 
function from each actuator to its associated sensor determines the relative level of 
control authority over each mode. Experiment 1 had shown poor LMED performance, 
especially with respect to low frequency modes. Thus, the AGS torquers were tuned 
to provide proportional-plus- derivative control of the cluster of modes near 1 Hz. 

The middle LMEDs (A and B) provide rate feedback control of the 4.5 Hz modes, and 
the lower LMEDs were tuned to the 8 Hz modes. The following control filter and gain 
matrices were obtained. The experimental results are discussed in the following 
section . 


V = 1 9 x 6 y X A X B X C X d1 

base rate LMED accelerations 
gyros 

A = diagonal [A A A A A. A A.] 
c & cx cy cz ca cb cc cd J 

(14 x 14) 


where 


A = A = A 
cx cy cz 


= I -0 " a 2Al 

L 1 --J 


. r.i27 

= E = 

_-29. 


.1270 .0018731 


ca 


cb 


r° a 2 B "| 

L 1 -ct 1B J 

= r° "“20"] 

L 1 -“ic-l 


A . = 

dc 


r° - a 2o 

L 1 '“ID- 


545 -.2026 


J 


a 1A = -1.49245 a 2A = 0.75575 


a lc = -1.39408 a 2B = 0.72654 


a 1(j = -0.72373 c* 2c = 0.71971 


a ld = -0 ' 51073 a 2d = 0,6928 


B = diagonal [B V B B B B B B ] 
c & cx cy cz cA cB cC cD 

(14 x 7) 
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where 


B = B = B - F = 
cx cy cz 


5.534E-5 
1. 8728E-3. 



"e lA 1 

"8.468E-5" 

B = 

in _ 


cA 

- 3 2A-I 

. 4. 687E- 3 . 


B 


cB 


3 IB 


" 9. 3962E- 5* 

- g 2B- 


. 5.178E-3 


B 


cC 


B 


cD 


” 3 ic" 


■-2.4993E-5' 

- 3 2C- 


_ 4. 4633E-3 . 

Pid' 


“ -3. 957E-4" 

Ud- 


_ 4. 496E- 3 


C c = diagonal [C cx C cy C oz C cA C cB C cC C oD ] 
(7 x 14) 

where 

C = [K G E] = [4.667E7 3.198E5] 

cx rx 

C cy = [K ry G E] = [5.599E7 3.838E5] 

C cz = [K rz G = [9. 332E6 6.396E4] 

C cA = [ ° K rA ] = [0 " 200] 

C cB = [0 K rB ] = [ ° - 200] 
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C eC = 10 K r C l = [0 - 901 

C cD = t» K rD > = - 100 ' 


DC = diagonal [D cx D cy D cB D( , c D cD ] 

= diagonal [ (Kpx * K^GF^ + K^GF^K^ + K^GFy- B oA ),<- S^)/- ^ )jC- 6 oD 

In the above definitions, matrices E, F, and G are those derived for a previous 
GF/LSSCV experiment; G = [1.58E4 1.58E5] while E and F are defined for matrices 
A c and B c . The gains K . represent rate feedback gains for the ith channel, and 

K . is the corresponding position gain; the following values were used: K^ x = 1000, 

K = 1200, K = 200, and K = -10, K = -12, K = -2. The remaining 
py pz rx i y rz 

parameters (o^, a^, 3 Q , f$2> K r ) f° r LMEDs A, B, C, and D were identified at 

the beginning of this section. Basic control equations are given in Table 4. 

IV. EXPERIMENTAL RESULTS 


This section presents selected portions of experimental results which are indica- 
tive of the performance obtained with the LMEDs. Results for each of the three con- 
trollers are presented separately. Section V contains some analysis and conclusions 
applicable to the data. The frequency response plots were obtained with an HP- 5423 
Dynamic Analyzer using bandlimited white noise for the input commands. Time history 
data was recorded with the HP-9000/COSMEC control computer and an analog strip- 
chart recorder. 


4.1 Experiment IB: MSFC Digital Controller 

Obtaining results during this experiment required more time than expected. 
Hardware problems in the control computer led to further delays. To interface the 
VCOSS hardware with the GF/LSSCV control computer, new software had to be written 
and verified, and the I/O capacity of the COSMEC had to be increased. These soft- 
ware changes were made along with a major reconfiguration of the main program, con- 
verting previously used assembly language segments from the COSMEC into HP-Basic. 
The additional I/O channels created intermittent problems with the current drive on 
one of the A/D converters. Meaningful data was nonetheless obtained as shown in 
Figure 15. The test configuration was with the AGS gimbals unlocked, and with the 
GF/LSSCV control computer active. Trial and error, with some simulation and analy- 
sis, led to the control gains and frequencies listed in Section 3.2. Again, very little 
performance was obtained from LMED-D. Figure 15a shows the open-loop response 
of each accelerometer to the excitation source at 4 Hz. Figure 15b shows performance 
gains with only LMED-A operating; notice the slight coupling into channel B. Figure 
15c shows the analogous result with only LMED-B operating closed- loop. Finally, 
Figure 15d shows the effect of operating all four LMEDs. Since the lower devices 
(C and D) were tuned to 12.5 Hz, they contributed little to this case. LMED-C 
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TABLE 4. BASIC CONTROL EQUATIONS 


Basic Control Equations 


Basic Equations: 


(2npxl ) 

Ap 

(2npx2np) 

Xp + 

(2npxl) 

Bp 

(2npxt) 

lip 

(*°1) 

Xs = 
(nsxl) 

A s 

(n s xn s ) 

X s + 
(nsxl) 

B S , 

(nsxm) 

Y P 

(mxl) 


(*x*) 

U„ + 

(AxS) 

Bu 

(Axt) 

(fel) 

(mx?) 

(mxinp) 

x p + 

(2npxl ) 

Dp 

(mxt) 

(Axl) 

Xc j = 
(ncxl) 

(ncxnc) 

XCi + 
(ncxl) 

Be 

(ncxns) 

Xs -j 

(nsxl) 


Discrete plant, sensor, actuators: 


wit! 


A Xi + B Ui 

(ndxl) (ndxl) (ndx*) (*xl) 


where nd * 2np + ns + * X-jT = [Xpj T X S i T Dp,!] 

2np ns l 


Xs-j = C Xi 

(nsxl) (nsxnd) (ndxl) 


+ 0 Vi ; C = [0 I ns 0] 

(nsxt)(txl) 2np nu 


Controller equations: 

Xci+1 

(ncxl) 


Ac Xci + Be Xs i = Ac Xc + Be 


Xi 


(ncxnc)(ncxl) (ncxns) (nsxl ) 


(nsxns)(nsxnd)(ndxl) 


where nc = 2ns = 2m 


Yci = Cc Xci + Dc Xsi 
(txl) (Axnc) (ncxl) (txns) (nsxl) 


Augmented state equations: 


li = C Xi 

(ns+ncxl) (ns+ncxnd+i+nc) (nd+A+ncxl) 


Hi = K Yi 

(Ixl) (txns+nc) (ns+nc) 


Xi T = C(X p i T Xs i ^ U Pi T) ViT Xc n -T] JjT = [Xs,^ XcJ] 

2np ns * *• nc (ns+nc) ns nc 
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Figure 15b. Response with LMED-A only closed loop. 










could effectively damp the 12.5 Hz mode, but in the process destabilized other modes. 
At this point, a linearity test was made using the LVDT on LMED-E (the device most 
visible from the control room). Figure 16 shows the open-loop response of LMED-E 
to low amplitude (0.1 V peak) triangular input signals at various frequencies below 
1 Hz. Figure 17 shows the response to constant frequency (0.1 Hz) varying amplitude 
triangular commands. Non-linear behavior is most obvious. At this point, the experi- 
ment was discontinued so that the devices could be analyzed more closely on a lab- 
oratory bench. Results from the device characterization tests are documented in a 
separate report. 


OPEN LOOP RESPONSE OP LMED E: CONSTANT AMPLITUDE (0.1 VPEAK) 

VARYING FREQUENCY (0.1 TO 1.0 Hz) 


AMP. 0.1 


0.1V PEAK 0.2 


0- 3 0.4 0.5 0.6 0.7 0.1 0.9 1.0 0.1 


LVDT RESPONSE 



INPUT SIGNAL 




Figure 16. LMED-E linearity check: LVDT response to constant 
amplitude, varying frequency. 


OPEN LOOP RESPONSE OF LMED-E: CONSTANT FREQUENCY (0.1 Hz) 

VARY AMPLITUDE (0.1 to 1.0 

AMP 0.12 0.14 0.16 0.1S 0.2 0.26 0.3 0.6 



Figure 17. LMED-E linearity check: LVDT response to constant 
frequency, varying amplitudes. 
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4.2 Experiment 2: MSFC I MED and AGS Controller 

Experiment 1 determined that the LMEDs were not operating in a linear fashion. 
At that point, the lower pair of devices (C and D) were removed from the test 
article and examined on a laboratory bench. In addition to mechanical misalignments 
and bearing contamination, the devices exhibited significant hysteresis and bearing 
roughness. At the same time, program schedules for future GF/LSSCV experiments 
required that the VCOSS experiments come to a conclusion. 

It was thus decided to use the spare LMED as an engineering model to reduce 
the dominant nonlinearities. The other LMEDs were used, as delivered, to demon- 
strate simultaneous operation of the AGS torquers and LMEDs in a decentralized con- 
trol scheme. Schedule constraints permitted only a day and a half to develop and 
demonstrate this complement of control hardware. Thus, a largely cut and try 
approach was undertaken. 

Figures 18 through 20 are sufficient to explain the procedure used. The con- 
trol algorithms were predetermined as discussed in Section 3.3; only the control gains 
and frequencies need be determined. Figure 18 shows the open-loop frequency 
response between LMED- A and its accelerometer. Clearly, the mode at 4.4 Hz should 
be controlled by the actuator. The previous results, however, demonstrated that this 
control effort destabilized other modes. For example, Figure 19 shows ihe relationship 
between the y-axis torque motor and accelerometer A; the upper response was 
obtained with LMED-A in the open-loop condition and the lower response with LMED-A 
operating closed-loop. The 4.4 Hz mode is virtually eliminated, but the 1 Hz mode 
has been destabilized. Not shown here is the destabilizing effect on a nearby 8.5 Hz 
mode. The y-axis torque motor was thus commanded to squelch the low frequency 
modes (below 2 Hz). The upper trace in Figure 20 shows an expanded open-loop view 
of that of Figure 19. The lower response incorporates both LMED-A and y-axis torquer 
closed-loops. Comparing Figures 19 and 20 the lower frequency modes have been 
actively damped, but the 8.4 Hz mode has been destabilized. The open-loop response 
at LMED-C would suggest sufficient control at 8 Hz to reduce this response as before. 
However, the device non-linearities prevented success in that approach. Figure 21 
shows the performance gains obtained by LMEDs A and B independently at their con- 
trol frequencies. Figure 22 shows the AGS control effectiveness at 1 Hz with the 
LMEDs also operating; damping is effectively doubled by visual inspection. Figure 
23 provides a dramatic demonstration of the effectiveness of LMED-B over a 5 Hz 
mode with all of the other actuators operating. (This response was obtained during 
a demonstration for the MSFC Workshop on Structural Dynamic and Control Interaction 
of Flexible Structures in April 1986.) Finally, the transient response of this control 
system is seen in Figure 24 by comparing open and closed loop responses to step 
commands on the Base Excitation Table. 


V. CONCLUSIONS 


1) The VCOSS- II LMEDs have the demonstrated capability to actively increase 
the damping of the LSS test article, using narrow band control filters. Experiment 
IB shows 60 percent increase in the damping of the controlled modes. In Experiment 
2 the LMED virtually eliminated the controlled vibration mode at 5 Hz. 

2) Non-linearities of the LMEDs make broadband control applications difficult. 
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0.0 HZ 12.000 

Figure 18. Magnitude of LMED-A to Acc-A transfer function in 
configuration 2; open loop. 
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Figure 19. Magnitude of torque-Y to Acc-A transfer function in 
configuration 2; (open loop) LMEDA. 

Closed loop effects of LMED-A on Y-torquer to Acc-A transfer function. 





Figure 20. Open and closed loop autospectrum of Acc-A with 
AGS and LMED-A controller. 
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Figure 21a. Closed loop response at 4.4 Hz; LMED-A controller only 
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Figure 21b. Open loop response at 4.4 Hz; LMED-A controller only 
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Figure 22a. Open loop response of rate gyro; BET at 1 Hz 
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Figure 22b. Open loop response of LMED accelerometer; BET at 1 Hz 
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Figure 22c. Closed loop response of rate gyro with all 
LMEDs and AGS; BET at 1 Hz. 
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Figure 22d. Closed loop response of LMED accelerometer with all 
LMEDs and AGS; BET at 1 Hz. 
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Figure 23a. Response of rate gyro at 5 Hz with all 
actuators except LMED-B. 
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Figure 23b. Response of LMED accelerometer at 5 Hz with all 
actuators except LMED-B. 
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Figure 23c. Response of rate gyro at 5 Hz with all actuators operating. 


1.0 

.8 

A .6 
.4 

B .2 
-.0 
C-.2 
-.4 
D -.6 
-.8 
- 1.0 

0.0 .3 .8 .9 ^2 1.5 1.8 2.1 2.4 2.7 3.0 



■ 


n 


n 

m 

a 

mi 


■ 

B 

Hi 

g 

■ 



■ 


m 

■ 

■ 

■ 

■ 




■ 


m 




m 

H 



B 


i 




m 

m 



B 


m 

m 


ilH 

m 

i 



B 


m 



m 


gg 





m 

B 

INI 

B 

3 

Nl 



■ 


m 


m 



KM 


! 


iH 

a 

M 

■ 

1 

m 

■ 



3 

a 


Figure 23d. Response of LMED accelerometer at 5 Hz with all 

actuators operating. 




Figure 24a. Open loop transient response of rate gyro (base excitation) 
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Figure 24b. Open loop transient response of LMED accelerometer 

(base excitation) . 
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Figure 24c. 


Closed loop transient response of rate gyro 


(base excitation). 
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Figure 24d. Closed loop transient response of LMED accelerometer 

(base excitation). 




3) The above comments require that a structure's open-loop frequency 
response be well known, either analytically or experimentally, prior to control law 
design with the present hardware. 

4) The LMED nonlinearities, and the attendant generation of frequency har- 
monics, complicate the analysis of experimental results. It is not clear whether the 
nonlinearities or the "control spillover" is responsible for the destabilizing tendencies 
of the LMEDs toward the "uncontrolled" modes. 

5) Although the supporting data is not presented, test data has shown that 
the middle LMEDs (A and C) and the AGS together achieved better vibration suppres- 
sion than either system could alone, using the given forms of control algorithms. 


VI. RECOMMENDATIONS 


1) The LMED hardware needs to be carefully reworked to eliminate the domi- 
nant non-linearities. Careful alignment of the linear shaft must be ensured. Further- 
more, the ball bearings should be replaced with high quality, preferably non- 
magnetic, components. The linear shaft should be replaced with high quality non- 
magnetic material. The reworked LMEDs should then be carefully characterized and 
supported with analytic models. See the report entitled "Characterization and Hard- 
ware Modifications of LMEDs" prepared by Control Dynamics for MSFC. 

2) Subsequent experiments should be planned carefully, allowing sufficient time 
to resolve unexpected problems which invariably occur. These experiments should 
include a period of thorough open-loop testing to validate analytic models. 

3) The control laws developed for future tests should first be designed, tested 
and evaluated based on the analytic models. These algorithms could subsequently be 
re-tuned based on experimental results. This would be more representative of the 
LSS control problem, and would address the requirements for on-orbit dynamic testing 
and verification. 

4) Detailed analytic predictions should be available to help analyze the experi- 
mental results. Such analysis might include eigenvalue and eigenvector analyses along 
with digital simulation capabilities. Such work was initiated at MSFC to support the 
VCOSS-II experiments. However, the results obtained only crudely agreed with the 
experimental data. This was attributed primarily to the unverified dynamic models 
and the unmodeled actuator nonlinearities. 
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